Repelling Random Walkers in a Diffusion-Coalescence System 



F. H. Jafarpour % and S. R. Masharian 2 
1 Physics Department, Bu-Ali Sina University, Ramadan, Iran and 
2 Institute for Advanced Studies in Basic Sciences, Zanjan, Iran 
(Dated: February 19, 2008) 

We have shown that the steady state probability distribution lunction of a diffusion-coalescence 
system on a one-dimensional lattice of length L with reflecting boundaries can be written in terms 
of a superposition of double shock structures which perform biased random walks on the lattice 
while repelling each other. The shocks can enter into the system and leave it from the boundaries. 
Depending on the microscopic reaction rates, the system is known to have two different phases. We 
have found that the mean distance between the shock positions is of order L in one phase while it 
is of order 1 in the other phase. 
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Recently the investigation of the microscopic structure 
and dynamics of shocks defined as discontinuities in the 
space dependence of the densities of particles in one- 
dimensional driven diffusive systems, has drawn much at- 
tention [lj]-[ll|]. It has been shown that the steady states 
of some of these systems can be explained in terms of col- 
lective excitations with one or more conservation laws. 
In [H three families of single-species driven-diffusive sys- 
tems are studied in which a traveling shock with a step- 
like density profile exists and behaves like a one-particle 
excitation in the system provided that the microscopic 
hopping rates are fine tuned. This has also been ob- 
served in the systems with more than one species of par- 
ticles @-[l(|. O n the other hand, the steady states of 
these systems can be written in terms of a superposition 
of such product shock measures. In @ and [ll| the au- 
thors have shown that such steady states are associated 
with the existence of two-dimensional representations of 
the quadratic algebras of these systems when they are 
studied using the Matrix Product Formalism (MPF) (for 
a recent review see |14j ) . According to this formalism the 
steady state of some of one-dimensional driven-diffusive 
systems can be written in terms of products of noncom- 
muting operators which satisfy a quadratic algebra. 
However, little is known about the microscopic dynamics 
of multiple shocks in these systems. The only example 
is given in [3| where multiple shocks are studied for the 
partially asymmetric simple exclusion process with open 
boundaries. In this paper we investigate the dynamics of 
a double shock structure in a branching coalescing sys- 
tem with nonconserving dynamics and reflecting bound- 
aries. The steady state properties of this system has 
already been studied in [Pa ] and (l6| . It turns out that 
depending on the microscopic reaction rates of the sys- 
tem it can be in two different phases: a high-density and 
a low-density phase. Since the dynamics of the system 
is non-conserving, the mean density of the particles in 
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the system in high-density phase is greater than that in 
the low-density phase. However, it has been shown that 
if one considers a canonical ensemble in which the total 
number of particles is conserved then the system has two 
phases: a high-density and a shock phase. In this case 
the shock does not have any dynamics [T^. In 0] the 
authors have shown that a single shock with biased ran- 
dom walk dynamics can evolve in this system provided 
that the boundaries are open so that the particles can 
enter and leave the system from there. Later in [l2| and 
[l3| it was shown that in an infinite system double shock 
structures with random walk dynamics can also evolve 
in the system. However, nothing is known about the dy- 
namics of these double shock structures in a system with 
boundaries. Our main attempt in this paper is to study 
the microscopic dynamics of such structures on a lattice 
of finite length and reflecting boundaries. 
In what follows we first define the model and then us- 
ing the Hamiltonian formalism show how a double shock 
product measure evolves in time under the Hamiltonian 
of the system. From there we construct the steady state 
probability distribution function of the system as a lin- 
ear combination of such double product shock measures. 
The mean distance between the shock positions is also 
calculated in the thermodynamic limit. 
The system in question consists of identical classical par- 
ticles on a one-dimensional lattice of length L. There is 
no injection or extraction of particles at the boundaries. 
The reaction rules between two consecutive sites k and 
k + 1 on the lattice are as follows: 
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in which A and stand for the presence of a particle 
and a hole respectively. As can be seen, the parameter 
q determines the asymmetry of the system. For q > 1 
(q < 1) the particles have a tendency to move in the 
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leftward (rightward) direction. For any q the model is 
also invariant under the following transformations: 
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in which k is a given site on the lattice. Throughout 
this paper we will only consider the case q > 1. The 
results for the case q < 1 can easily be obtained using 
(|2|). By formulating the stochastic Hamiltonian of the 
system as a quantum spin chain, it has been shown that 
it is completely integrable [HI, [ly] . As we mentioned the 
system has two different phases depending on the values 
of q and A. For q > 1 it has a high-density phase for 
q 2 < 1 + A and a low-density phase for q 2 > 1 + A. 
In the high-density phase the density profile of particles 
has its maximum value near the left boundary while it 
is a constant p = in the bulk of the lattice. It 

also drops exponentially to zero near the right boundary. 
The particle correlations exist at both boundaries. In the 
low-density phase the density profile of particle has again 
its maximum value near the left boundary but it quickly 
drops exponentially to zero in the bulk and remains zero 
throughout the lattice. In this phase the particle cor- 
relations only exist near the left boundary. The mean 
density of particles is of order j- in this phase. On the 
transition line q 2 = 1 + A the density profile of particles 
drops exponentially near the left boundary while changes 
linearly in the bulk of the system. The mean density of 
particles in the bulk of the lattice is equal to 2 (i+a) m 
the thermodynamic limit. 

Recently, it has been shown that the steady state prob- 
ability distribution function of some of one-dimensional 
driven-diffusive systems can be written in terms of inter- 
actions of single shock structures @ . In the following we 
will show that the steady state of our coalescence system 
defined by ([T]) can also be written in terms of superposi- 
tion of double shock structures. These shocks repel each 
other while perform biased random walk on the lattice. 
Any state of the system is defined through a proba- 
bility measure P^ on the set of all configurations r\ = 
(Vi,il2] • ■ ■ i Vl), Vk £ {0; A}. F° r our purposes it is con- 
venient to use the Hamiltonian formalism where one as- 
signs a basis vector \rj) of the vector space (C 2 )® L to each 
configuration and the probability vector is defined by 
\P) = -^l 7 ?) which is normalized such that (s\P) = 1 
where (s\ = XL^I- The time evolution is now described 
by the master equation: 
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in which H is called the Hamiltonian and its matrix ele- 
ments are the hopping rates between any two configura- 



tions. For a system defined on a lattice of length L with 
reflecting boundaries the Hamiltonian can be written as: 
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FIG. 1: Sketch of a double shock structure. The shock posi- 
tions are defined at the sites m and n.. 



where h^^+i acts nontrivially only on sites k and k 
In a basis defined as: 
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the local Hamiltonian of our system in (dJ has the fol- 
lowing form: 
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We define a double Bernoulli shock measure which is a 
product measure with two jumps in the local particle 
density associated with two random walkers (the shock 
fronts) at sites m and n as: 
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in which < m < n — 1 and 1 < n < L + 1. Here we 
have introduced two auxiliary sites and L + 1. A sim- 
ple sketch of such shock measure is given in FIG [TJ It is 
easy to verify that this family of shock measures generate 
a subspace of the vector-space of states which is invari- 
ant under the time evolution generated by H and thus 
the many-particle problem is reduced to a two-particle 
one. As we mentioned earlier, the time evolution of such 
product shock measure has already been studied for an 
infinite system with no boundaries [12, [HI; nevertheless, 
in this paper we aim to study a finite system with reflect- 
ing boundaries. The time evolution equations for \P m ,n) 
are given by: 
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H\P m ,n) = q- x \P m+ x,n) + 9(1 + A)|P m _i )B > + q-^l + A)|P m , n+1 ) + q\P m , n - X ) - (q + q~ x ){2 + A)|P m , n ) for 
m = 1, • ■ • ,L — 2 and n = m + 2, ■ • ■ , L 

P|Po.„) = (T 1 - <z)|Pi.„) + + A)|P ,„+i> + 9|P ,n-i> - ? -1 (2 + A)|P ,„) for n = 2, • • • ,L 
H\P m ,L+i) = g -1 |Pm+i,L+i> + 9(1 + A)|P ro _i, L+ i) + (g - g-^lPm.z) - ?(2 + A)|P m x+i> for m = 1, 1 
H\P 0>L+1 ) = (q- 1 - q)\P hL+ i) + (q- q' 1 )^) 
#|P m , m +i) = for m = 0, 

I 
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As can be seen for q > 1 the left random walker per- 
forms a biased random walk and preferentially hops to 
the left regardless of the values of q and A. In contrast 
the right random walker preferentially hops to the left for 
q 2 > 1 + A and to the right for g 2 < 1 + A. On the co- 
existence line where q 2 = 1 + A the right random walker 
performs an unbiased random walk. The left (right) ran- 
dom walker can also leave the lattice only from the left 
(right) boundary. The diffusion coefficients and also the 
velocities of the random walkers can now be easily calcu- 
lated from ©. 

Let us now explain why the random walkers repel each 
other. It can easily be seen from JSJ that as long as the 
shock positions are more than a single site apart, they 
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never meet each other during the time evolution. How- 
ever, it seems from there that the random walkers can 
meet each other when they are a single site apart. In 
what follows we show that this is not the case. For in- 
stance we consider the first equation in JSJ for n = m + 2 
where the shock positions are a single site apart. Rewrit- 
ing this equation in terms of a new definition for the 
shock measure as: 
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As can be seen the shock positions never get closer that 
a singe site. In fact the dynamical rules ([T]) do not allow 
the shock fronts to get closer than a single site since it 
results in an empty lattice. This is why we say that the 
random walkers repel each other. One can easily check 
this for other equations in ([8]) in which the shock posi- 
tions are a single site apart to see that in terms of the 
\P m ,n} the shock positions never meet and the minimum 
distance between them is at least a single site. 
In this paper we are specially interested in the steady 
state of the system. One should note that an empty 
lattice is a trivial steady state for the system. It can be 
seen from ([IJ that an empty lattice never evolves in time. 
There is actually a nontrivial steady state for the system 
in which the lattice contains some particles. The non- 
trivial steady state of the system can now be constructed 
as a superposition of double shock measures as follows: 
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provided that we exclude the empty lattice from \P*) by 
requiring: 
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in which: 



|0> = |0} 0L = 
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is associated with a configuration with no particles in the 
system. The normalization factor Zl in ITU)) can easily be 

obtained from Z L = J2 m =o E n =l+i i>m,n- By requiring 
H\P*) — we find a system of equations for V'm.n's. It 
turns out that this system of equations has the following 
solution: 
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Note that ipm,n's m (fl3|) are also valid for g 2 = 1 + A and 
at the coexistence line one should only replace A with 
q 2 — 1 in (flU|) . One can see from ([7]) that there are L + l 
states in which the shock positions are at two consecutive 
sites. The states |P mim +i)'s point to an empty lattice. 
Since the empty lattice is a trivial steady state of the 
system; therefore, the coefficient of these states in (fTU|) 
i.e. "i/vn.m+i's, are taken to be equal to ip' . The condition 
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It turns out that on the transition line q 2 
becomes: 
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The normalization factor which is called the grand- 
canonical partition function of the system can now be 
calculated and after substituting ip' from (| 14[) and (fT5l) 
one finds: 
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for <? 2 = 1 + A. 
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As one can see our results obtained here are exactly those 
obtained in [lj| and [l6| using different approaches. Us- 
ing the steady state probability distribution function (flQ|) 
one can easily calculate the density profile of the particles 
and also any correlations in the steady state. However, 
since the results are exactly those obtained in the above 
mentioned papers, the results are not given here. 
Having the probability of finding the random walkers at 
sites m and n in the steady state, one can calculate the 
mean distance of the shock fronts in the steady state de- 
fined as: 
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It turns out that in the thermodynamic limit L 
is given by: 
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In the high-density phase q 2 < 1 + A the shock fronts 
have their maximum distance while in the low-density 



phase they have the minimum distance which is of order 
of a single site. One should note that the mean distance 
of the two shock fronts changes abruptly from one phase 
to the other phase which can be a sign for the phase 
transition in the system. 

It is also interesting to study the probability of finding 
each shock front at a given site in the steady state. The 
probability of finding the left shock front at the site m is 
defined as: 



L + l 



Pm = E 



for < m < L. 



(19) 



n— m+1 



In the thermodynamic limit L — * oo and in the high- 
density phase P m is an exponential function with the 
inverse length scale ln(g 4 ) while in the low-density phase 
it is an exponential function with the inverse length scale 
equal to ln(q 2 (l + A)). On the other hand the probability 
of finding the right shock front at site n is given by: 



P n = ^2 V-Vn for 1 < n < L + 1. 



(20) 
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In the thermodynamic limit L — > oo and in both the 
high-density and the low-density phase this probability 
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distribution function has an exponential behavior with 
the inverse length scale ln(g 2 (l + A) -1 ). This explains 
why the system has three diffrent length scales. 
In this paper we have studied a coalescence system with 
reflecting boundaries and showed that its steady state 
can be explained in terms of superposition of probabil- 
ity distribution of two interacting random walkers which 
perform biased random walks while repelling each other. 
The random walkers can also leave or enter from the 
boundaries. One should note that the random walk pic- 
ture actually fails at the left boundary for q > 1 (and at 
the right boundary for q < 1). In fact, as can be seen 
from ||5J), the left random walker should enter the system 



with a negative rate. This has already been observed 
in the branching-coalescing model with open boundaries 
studied in [4|. Apart from this, we have found that the 
steady state probability distribution function of the sys- 
tem is exactly the one obtained in (l5l . [l6j which ob- 
viously generates the same density profile of particles in 
the system in each phase as it was calculated by the same 
authors. It is interesting to consider a more general re- 
action rates in JT]) and see under what constraints the 
random walk picture in a system with refelecting bound- 
aries exists. This is under our investigations and will be 
published elsewhere. 
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